Entanglement criterion for pure M N bipartite quantum states 
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We propose a entanglement measure for pure M ® N bipartite quantum states. We obtain the 
measure by generalizing the equivalent measure for a 2® 2 system, via a 2® 3 system, to the general 
bipartite case. The measure emphasizes the role Bell states have, both for forming the measure, and 
for experimentally measuring the entanglement. The form of the measure is similar to generalized 
concurrence. In the case of 2 <g) 3 systems, we prove that our measure, that is directly measurable, 
equals the concurrence. It is also shown that in order to measure the entanglement, it is sufficient 
to measure the projections of the state onto a maximum of M(M — 1)N(N — l)/2 Bell states. 

PACS numbers: 03.67.Mn,42.50.Dv,42.50.Hz,42.65.Ky 



I. INTRODUCTION 

The concept of quantum entanglement is not new, it 
goes back to the early days of quantum theory where it 
was initiated by Einstein, Podolsky, and Rosen Q], and 
Schrodinger 0. Many years has passed since the dawn 
of quantum mechanics, but we have still not been able 
to solve the enigma of entanglement, e.g., finding a com- 
plete mathematical model to describe and quantify this 
interesting feature of quantum mechanical systems, and 
in the same time reveal the physical implications of this 
feature. During recent years, separability and entangle- 
ment has been a vital research field. Peres pioneered 
quantification of entanglement by showing that a neces- 
sary criterion for separability was positivity of the density 
matrix upon partial transposition 0. Soon thereafter, 
Horodeckis proved that the criterion was also sufficient 
i. 

Recently, several quantitative measures of entangle- 
ment for nonseparable states, such as entanglement of 
formation [5(, distillable entanglement |fj, relative en- 
tropy of entanglement H, concurrence , or concurrence 
related measures 0, llflllll Il2l IT^j have been suggested. 
In particular, the definition of concurrence is based on 
the spin-flip operation, and Rungta et al. ^3 have gen- 
eralized this operator and defined I-concurrence for pure 
bipartite state in any dimension. For the mixed 2(g) M bi- 
partite state we have a separability criterion given in [l4| 
and a lower bound of concurrence of mixed such quantum 
states 0| • Entanglement witnesses is another method of 
detecting entanglement 

In this paper, we develop a measure for general pure 
N®M bipartite quantum states, inspired by the measure 
of entanglement proposed in |17lll8l |. It is based on bipar- 
tite phase sums and differences, and we conjecture that if 
it is properly normalized, it equals the I-concurrence for 
all pure bipartite states [U$. In Sec. ITTlwe briefly discuss 



the physical grounds for our measure for 2 <g> 2 systems 
(two qubits). In Sec. IHII we argue what a similar measure 
should be for a 2 ® 3. In Sec. IIVI we derive a measure 
using a relative-phase positive operator valued measure 
(POVM) as the entanglement quantifier. In Sec. Ivl we 
compare our measurement to existing measures. In Sec. 
I VII we extent the measure to encompass bipartite systems 
of any dimension. In Sec. I VIII we apply our measure to a 
few sample states and discuss its normalization and the 
role of the Bell states. Finally, in Sec. IVIIII we summa- 
rize our findings. The main novelty with our paper is not 
the measure itself, since it is essentially proportional to 
I-concurrence, but rather the way the measure is derived. 
We have based our derivation on physical argumentation 
in contrast to earlier derivations that are primarily based 
on mathematical arguments. 



II. ENTANGLEMENTBETWEEN TWO 
QUBITS, A BRIEF REVIEW 

In an earlier paper we investigated the entanglement 
properties of a 2 ® 2 bipartite state 0| . The starting 
point of that investigation was the assumption that the 
entanglement properties are, or can be, expressed in the 
bipartite state's joint phase properties. On basis of this 
assumption, we found that the state's non-local prop- 
erties are found in two of the state's off-diagonal den- 
sity matrix coefficients (when the density matrix was ex- 
pressed in the standard basis |11), |12), |21), and |22)). 
Examining the reasons for this, a rather simple physical 
reasoning shows why this is the case. Tensor multiply- 
ing two qubit density matrices, p = pa <£> Pb, where each 
factor pi, i = A, B, is of the form 
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and subdividing the ensuing 4x4 matrix into four 2x2 
quadrants, it is clear that the upper left quadrant does 
not contain any information about qubit A's phase ifA- 
The lower right quadrant lacks this information, too. Due 



2 



to hermiticity, the upper right and the lower left quad- 
rant contain the same information. Therefore it suffices 
to consider, e.g., the upper right quadrant. Of its four 
coefficients, its diagonal terms (that is, the coefficients 13 
and 24 of the joint system density matrix) do not con- 
tain information about qubit B's phase ipB- Hence, all 
the joint phase information is collected in the off-diagonal 
coefficients 14 and 23, and we found that a relevant mea- 
sure of entanglement for a two qubit system was 



Tsup = Sup(2||pu| - |/323| |), 



(1) 



where Sup refers to the supremum of the function with 
respect to any local unitary transformation(s). 

One then notes that by a local phase rotations, it is 
always possible to make ,014 and P23 simultaneously real. 
In this case 

Re(pi 4 ) = |/314| = (*+|/3|*+> - (*-|p|*->, (2) 



Re(p 2 3) - |P2 3 | = (*+|p|*+) - (*-|p|*- 



(3) 



where |*±) = (|11) ± |22))/V2, |$±) = (|12) ± |21))/V2, 
and Re denotes the real part. That is, the entangle- 
ment is simply the maximum of the difference between 
the state's projections onto the Bell states. Hence, it can 
be measured as the Bell-analyzer visibility [Tsl ]. 



III. AN ENTANGLEMENT MEASURE FOR A 
2x3 BIPARTITE STATE 

Now we apply a physical reasoning similar to that in 
Sec. [H]to a bipartite 2 <g> 3 system (a qubit and a qutrit). 
We start by subdividing the system density operator (a 
6x6 matrix) into four 3x3 matrix quadrants. It is 
clear that the upper left and lower right quadrants do 
not contain information about the qubit phase </ja- The 
remaining two quadrants contain the same information 
so let us focus on, e.g., the upper right quadrant. Disre- 
garding the diagonal coefficients in this quadrant (they 
contain information about ip& but not about the qutrit 
phases), we see that the joint phase properties of the 
state lies in the density matrix coefficients 15, 16, 24, 
26, 34, and 35. These coefficients correspond to the 
projectors |11)(22|, |11)(23|, |12)(21|, |12)(23|, |13)(21|, 
and |13)(22|, respectively. Consider now the two possi- 
ble complete Bell bases: 



l*i> 
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(|11> + |22)), |* 2 ) 
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* 3 ) = -=(|12) + |23», |*4) = -=(|12) - |23» 



V2 



V2 V 



and 



|*5> - -=(|13) + |21)), |* 6 ) = -=(113) - |21)) 



V2 



V2 V 
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(|H) + |23», |* a 



V2 



(|11> - |23», 



|* 3 ) = ^(|12) + |21)), |$ 4 ) = ^=(|12)-|21)), 



and 

l$5> - ^(|13) + |22», |* 6 ) = -J=(|13) - |22)). 

The two bases can be obtained from each other by per- 
mutation of any two of the qutrit states. However, such 
a permutation changes all six states, that is, one whole 
basis is transformed into the other. We note that 

Pis +P51 = 2Re(p 15 ) = (*i|/5|*i) - (* 2 |/5|*2>, (4) 



and 



Pie + Pm = 2Re(p 16 ) = (*i|/3|*i) + (* 2 |/5|*2>, (5) 

etc. By applying phase shifts, local to the qutrit B, it 
is always possible to make all three coefficients P15, P26) 
and P34 real, simultaneously. Hence, the absolute values 
of density matrix coefficients 15, 26, and 34 are associ- 
ated with, and can be obtained from, a Bell-state analysis 
using the complete basis set {l 1 ^)}, while, in a similar 
manner, coefficients 16, 24, and 35 are associated with 
the complete and noncompatible set {|$i)}. Intuitively, 
one would expect the entanglement to be greatest if only 
the states in one of the sets were excited. Moreover, 
the entanglement should be maximized if only one of the 
states in one of the sets were excited. A reasonable mea- 
sure of entanglement would therefore be 



r = (M(||pib 

+ I|P34| - \P16 



P24 
2 



IP26 1 



|P35|f 



(6) 



where, again, T is not invariant to transformations local 
to qubit A and qutrit B, respectively, so the state's en- 
tanglement is understood to be given by Tsup, the supre- 
mum of r taken over all possible local unitary transfor- 
mations. A/2 is a normalization factor that, if taken to 
be 2, makes < T < 1. Eq. © is our central result for 
2 Cg) 3 systems. Due to the absolute signs, the expression 
is symmetric with respect to the two Bell basis sets. It is 
also clear that for any separable state, p = pa <£> Pb , one 
gets 



r = 



(7) 



Since P15 = PA12PB12, P24 = PA12/3 B 12' P26 = PA12/3B23, 
P35 = PA12PB23' an d PW = PA12/3B13, P34 = PA12Pb13' 

These relations explain why we have chosen the specific 
pairing between the coefficients in sets {P15, P2Q, P34} and 
{pi6, P24, P35} associated with the respective Bell basis. 
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IV. A MATHEMATICAL DERIVATION 

In order to eventually generalize the results, we need 
to be on a little bit firmer mathematical ground. We 
start by introducing the basis vectors |1)a, |2}aj and 
|1)b) |2)b, |3)b for subsystems A and B. Subsequently we 
form the Hermitian operator 

Aa^a ; i 2 ) = -^(/ + e^ 12 |l) A (2| A + h.c), (8) 

ZTT 

where h.c. denotes the hermitian conjugate. In the same 
manner we define 

A B (^B;12,VB;13,</>B;23) = ^ V + | 1) B (2 1 B 

Z7T 

+e^ 13 |l} B <3| B (9) 
+e^ B ^|2) B <3|B+h.c.). 

The bipartite system's phase properties are described by 
the operator 

A(<^A;12, ¥>B;12, V?B;13 , </?B;23 ) = A A ® A B . (10) 

We can re-express this operator in terms of the system's 
sum and difference phases 

VplU = ^A;kl ± fB-pq- (H) 

The linear dependence between i^A;kl> V?B;pq and V p ) g +, 
tPp'q- allows us to express the operator in (fTUf> as a func- 
tion of the sum and difference phases. Hence we can 

A 12 12 12 12 12 12 

write AO^, ip x \ 2 _, <p x ',3+, f {,3-^2,3+ ^2,3-)- Next we 
compute to what extent the density operator depends on 
these phase sums and differences. Since the sum and dif- 
ference phase POVM is a periodic function of the phases, 
we can compute the Fourier components of the POVM's 
expectation value. We define, e.g., 

L Jin 

rrw-r/W 1-2 1,2 1,2 1.2 1,2 1.2 \ 

+ Mv\' 2 2+ + 1") <Pi%- + A", <p{% + + 7T, tp\l_ + TT 
>¥>2\3+ +^, ( fil' 2 3- +7r)])|- (12) 

The addition of all the it terms in the right-most term in- 
side the trace operation above makes the functions T p ' q + 
above insensitive to the diagonal coefficients of the den- 
sity matrix, so only the joint (nonseparable) properties 
are probed by T p ' q ±. Finally, we compute 

r = (27rAA 2 [||r^ + |-|rl ; Lll 2 + l|r 1 ;^|-|r^_|| 2 
+l|r^ + l-|r 2 :Lll 2 ]) 1/2 

= (^[Hpisl-l^ll' + ll^el-lpssll 2 (13) 

+ IIP34|"|/5i 6 || 2 ]) 1/2 . 

The result is identical to the one we obtained by reason- 
ing in Eq. ((HJ), above. The derivation may seem rather 
ad hoc, but follows our derivation for the similar entan- 
glement measure for 2(8)2 systems. The various steps are 
motivated in [Tsl ]. 



V. A COMPARISON TO EXISTING MEASURES 

In this section we compare our degree of entanglement 
with concurrence for pure bipartite quantum systems. 
There are several attempts to generalize concurrence for 
pure bipartite quantum systems. One of these general- 
izations is called I-concurrence and is defined in terms of 
a super operator that is a generalization of a spin 1 /2 flip- 
operation to higher Hilbert-space dimensions. For a pure 
state |«&), of dimension TV eg) M, I-concurrence is defined 
[13 as 

d = ^2(1-Tr(e|)) = ^2(1-Tr(p 2 )), (14) 

where g A = Tr B (|$)($|) and g B = Tr A (|$)($|). Let |*) 
be a pure state for a 2 3 bipartite quantum systems 
given by 

23 23 

w = EE ® i*> = EE a ^ kl ) ( 15 ) 
fc=i 1=1 <s=i 1=1 

where {\k)} and {\l}} are two complete orthonormal ba- 
sis vector sets spanning a M- and iV-dimensional Hilbcrt 
space, respectively, and where J2k=i Sj=i \ a ki\ 2 — 1- 
For such a state, the I-concurrence Ci coincides with the 
concurrence C for the state an d it can be rather 
simply expressed in the state's probability amplitudes as 

C = (2(|aiia 2 2 - ai2a2i| 2 + |aii"23 - ai3«2i| 2 

+ |ai 2 a 23 - ai 3 a 22 | 2 )) 5 . (16) 

(Very similar expressions have been also derived in [TT| 
and in 01)- This expression should be compared to our 
entanglement measure 

Tsup = Sup[A/" 2 (|||aiia 22 | - |ai 2 a 21 || 2 

+ l|aii«23l - I"i3a2il| 2 

+ ||a 12 ^ 3 |-K 3 ^ 2 || 2 |)]i (17) 

It is seen that for A/2 = 2, we get Tgup < C. Suppose we 
rotate this state with local unitary operations Ua ® Ub 
such that 

2 3 

u A ®u B \i>) = J2J2 a 'u\ kl )- ( 18 ) 

fe=l 1=1 

The state's concurrence is invariant with respect of the 
operation. Hence, if it is always possible to make, e.g. 
a'n = a' 23 = 0, then the two expressions above coin- 
cide to yield (2(\a' 12 a' 21 \ 2 + \a' 13 a' 21 \ 2 + K 3 a' 22 | 2 ))5, and 
such a transformation maximizes T. We shall now prove 
that this is always possible. Consider the local unitary 
transformations 

U ^ 9 '^ = [iwxe encase ) 
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and 

I 1 ° ° \ 

U B (0,<p) = cos i? ie^sint? . 

\ isintf e^cos-d ) 

From these transformations we get 

a' n — an cosO + ia 2 ie % ^ sin#. (19) 

It is seen that it is always possible to make a' n =0 by 
an appropriate choice of 9 and <j>. We also get 

a 23 = * sini9(iai2 sin 6* + 0(226 cos 8)) 

cosi9(iai 3 sin (9 + a 23 e^ cos 9). (20) 

For a given choice of 9 and </>, the expressions in the 
parenthesis on the right hand side of the equation above 
are two fixed complex numbers Therefore, it is always 
simultaneously possible to assure that a' 23 = by an 
appropriate choice of the parameters $ and ip. Hence, the 
proof that our measure coincides with the concurrence for 
pure states is completed. 

From (|17|l above, we see that it is formally also possi- 
ble to prove that the measures coincide if, e.g., a' n — 

"23 = °> or if a ll = a 22 = °> a ll = «'l2 = °) or 

a'n = a' 13 = 0, where we have looked at the possi- 
bilities where a' n = 0. (The other possibilities lead 
to equivalent conclusions.) The first case was already 
proven, above. In the second proof outlined 

like the one pertaining to the first case, shown above, 
can be used. In the third case we can prove the asser- 
tion as follows: The pure state can be written \ip) — 
|1) <g) (au|l) + «i2|2) + ai3|3» + . . .. We see that in 
order to prove our assertions, it only makes sense to con- 
sider rotations in subspace B. However, in order to make 
a'n = a' 12 = we need to find three new basis vectors 
in the subspace that are mutually orthogonal, and that 
makes aii|l) +ai2|2) +ai3|3) — > a' 13 |3'). It is quite obvi- 
ous that the choice |3') = (an\l) + a u \2) + ai 3 \3)) / VJV , 
where Af is a normalization factor, and where |1') and |2') 
can be chosen arbitrary, as long as they are orthogonal to 
|3') and to each other, satisfies our requirement. In the 
fourth case we have \tp) = (an|l) + 0:21 1 2)) ® |1) + — 
We see that in order to get a' n = a' 21 = 0, we need 
to find a unitary transformation in space A rendering 
a n |l) + a 2 i|2) -> Oil') + 0|2'). In is obvious that if ei- 
ther an or «2i are nonzero, no such transformation can 
be found. 

As a consequence of the conclusions above, we can de- 
duce that in order to measure the entanglement of a pure 
state, one can do Bell measurements, but it is neces- 
sary to use not one, but a minimum of two different Bell 
bases. Suppose, e.g., that we have found a basis such 
that a' n = a' 23 = 0. We then need to find \p' 24 \, \p' 35 \, 
and |p 34 | to get r SU p. The first of the density matrix 
coefficients can be found by projection on the states |$i) 
and I $ 2 ) subtracting one of the outcome probabilities 
from the other, taking the absolute of the difference, and 



dividing by two. To obtain p' 35l we can likewise project 
on the states ($5) and (orthogonal to the first pair 
of projectors). Finally, to get p' 34 , we need to project on 
the states |^s) and l^e)- These two states, however, are 
not orthogonal to the any of the previous four Bell states. 

VI. GENERALIZING TO M <g> N BIPARTITE 
QUANTUM STATES 

The derivation of the entanglement measure for a 2 ® 3 
system made in Sec. IIVI can be extended to an ar- 
bitrary bipartite system. We label the basis vectors 
|1) A , \M) A , and |1) B , . • . , |iV) B . Subsequently we 
form the operator 

1 M-l M 

Aa = ^( /+ E E e^|fe) A (/|A + h.c), (21) 

k=l l=k+l 

where h.c. denotes the hermitian conjugate. In the same 
manner we define A B . Then we form the tensor product 
operator 

A = A A «A B . (22) 

We re-express this operator in terms of the system's sum 
and difference phases, given in Eq. (|llfl above, as 

Afal^+J^l^-J ■ ■ •jVjv_i,jv+' Vn-1,N- ( 23 ) 

1,3 1,3 M-l.M M-l.M \ 

1 ^l,2+> • • • ! 'Pjf-ljN—i ' ' ' ' ^1,2+ i • • • J i Pn-1,N-)' 

This operator depends on M(M - 1)N(N - l)/2 phase 
sums or differences. We define, e.g., 

= \\l dptt+j*'"* ( 24 ) 

rri / »r A I 1,2 1.2 M-l.M \ 

Tr(p[A(v? 1 ; 2+ , ip{ i2 _, <p N _x iN J 

a / 1,2 . 1,2 . M—l,M 1 mm 

+ A Wl]2+ + ^ <Pl,2- + 7r > • • " > V N -l,N- + n )\)\- 

Our final result, general to any bipartite system, is 
T = (2tt (25) 

M-l M N-l N 

^EEEE ll^ + |-|r^_|| 2 )^ 

fc=l l=k+l p—1 q=p+l 
M-l M N-l N 

= (^E EE E \\P(k-l)N+p,(l-l)N+q\ 
k=l l=k+l p=l q=p+l 

-|/5(fe-i)Af+,j,(;-i)Af+p|| 2 ) 1/2 - (26) 

This is our main result, where the entanglement Tsup of 
the state p is understood to be the supremum of the equa- 
tion (|26fl under all local unitary transformations. For a 
separable state p = px ® Pb, the expression (|2*S|l . above, 
simplifies to 

M-l M N-l N 

r = ^EEEE ( 27 ) 

k=l l=k+l p=l q=p+l 

||/5A;feiP B ; Pg | - \pA-kiP^;qp\\ 2 ) 1/2 = 0- 
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That is, our measure is identically zero for any separable 
state. Hence, the expression can be used as a separability 
criteria for any state, not only pure states. For a pure 
state where 



where K is the smaller of M and N. For this state we 



M N 



(28) 



fc=i 1=1 



we have 



1/2 



(Af-1 A/ AA-1 iV \ 
E I! I! E IKp<J - |afega* p || 2 J 
fc=l i=fe+l p=l g=p+l / 

(29) 

We can compare this expression with the concurrence (or 
similar measures such as the concurrence vector) for pure 
M <g> N bipartite states [1 El El 03 



Coc 



(M-l M N-l N \ 
E E I! I! \a kp ai q - a kq ai p \ 2 \ 
fe=l l=k+l p=l q-p+l j 



1/2 



(30) 

We conjecture that, properly normalized, Tsup = C for 
any state. An extensive numerical testing has always 
confirmed the hypothesis. We are presently working on 
a strict mathematical proof of the conjecture. 



VII. BIPARTITE ENTANGLEMENT, THE BELL 
STATES, AND NORMALIZATION 



get r 



sup 



a: 



1/2 



When M = N = K, such states 



are, e.g., manifested as the eigenstates of the relative- 
phase operato r Il9ll . and they have been experimentally 
demonstrated l20|. The states also have a role in quan- 
tum polarimetrv |2l| , and in quantum cryptography |22| . 
In spite of this, they seen to have no direct role as quan- 
tifying states of entanglement. 

The observation above raises an important question: 
How do we normalize Ts U p (or equivalently, the con- 
currence)? If we set A/2 = 1 we get < r SU p < 1. 
However, a Bell state in this Hilbert space will then 
have Tgup = (1/2) 1 / 2 < 1. Another possibility is to 
set A2 = 2 rendering Tsup = 1 for a Bell state and 
< Tsup < (2) 1 / 2 in general. Ideally, one would like 
to find a way to quantify entanglement so that one gets 
a quantitative measure that corresponds to the state's 
utility in quantum information tasks. Such a general 
quantification hierarchy lies beyond the scope of this pa- 
per. Investigations of minimum reversible entanglement 
generating sets 0, 0] have shown that even for rela- 
tively simple Hilbert spaces, such as those housing tripar- 
tite qubit systems, there exist two classes of states, W- 
states and GHZ-states, that cannot be transformed into 
the other class by local operations HSHH^}. Hence, it 
is not clear to us how an entanglement hierarchy should 
be defined. 



From the previous section we conjecture that by mea- 
suring certain of the density matrix coefficients, a pure 
state's entanglement can be quantified. All the needed 
coefficients can be obtained by projections on Bell states. 
E.g., the coefficient p^-i)N+p,(i-i)N+q can be obtained 
by projection onto the states (\kp) ± \lq))/y/2. We 
note that for an M <g> N bipartite system, there exist 
M(M - l)N(N - 1) Bell states. However, not all of 
those states are needed to measure the bipartite state's 
entanglement. In analogy with our results for 2 (g> 3 sys- 
tems, we conjecture that at most M (M-l)N(N- 1) /2 of 
the Bell states are needed, due to the degrees of freedom 
local unitary transformations give us. For M, N > 2 this 
implies that more than one Bell basis is needed, in gen- 
eral, to measure the state's entanglement (we can note 
that when the product MN is odd, complete Bell bases 
do not exist). At any rate, our results demonstrate that 
for bipartite states, the Bell states play a fundamental 
role in defining entanglement properties. 

At first sight, the central role of the Bell states in en- 
tanglement classification may seem obvious. However, 
perhaps surprisingly, these states do not have the largest 
entanglement of the bipartite states. For a Bell state, 
Tsup = (A2/2) 1 / 2 . However, the upper limit for Tsup is 
reached for, e.g., the state 

1 K 

— Y\kk), (31) 



VIII. DISCUSSION AND CONCLUSION 

In conclusion we have derived a quantitative measure 
of entanglement for pure M®N bipartite quantum states, 
based, essentially, on simple physical considerations. The 
measure is positive, bounded, invariant to local unitary 
operations, and it is shown that it equals zero for all 
separable states. Our measure is less than, or equal to 
generalized concurrence for pure bipartite state. We con- 
jecture that it is actually always proportional to the I- 
concurrence. In contrast to the latter, our measure can 
be obtained by measurement in a direct fashion. To this 
end, one projects the state onto Bell states and maxi- 
mizes the the sum of certain differences with respect to 
local unitary transformations. The fact only projections 
onto Bell states are needed suggests that these states are 
of particular significance for all bipartite systems. 
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